Non-Abelian Energy Loss at Finite Opacity 
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A systematic expansion in opacity, L/X, is used to clarify the non-linear behavior of induced 
gluon radiation in quark-gluon plasmas. The inclusive differential gluon distribution is calculated 
up to second order in opacity and compared to the zeroth order (factorization) limit. The opac- 
ity expansion makes it possible to take finite kinematic constraints into account that suppress jet 
quenching in nuclear collisions below RH1C (y/s = 200 AGeV) energies. 
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Introduction. The production of high transverse mo- 
mentum jets in QCD is always accompanied by gluon 
showers. For jets produced inside nuclei via e + A or 
quark-gluon plasmas via A + A, final state interactions 
of jet and radiated gluons induce further radiation and 
also broaden the gluon shower. The non-abelian radiative 
energy loss in a medium is expected to be observable as 
"jet quenching" . The observable consequences in nuclear 
collisions should be seen as a suppression of the high pj_ 
tails of single hadron distributions |l],^] and a broadening 
of the jet cone j|. 

Non-abelian energy loss in pQCD can be calculated an- 
alytically in two limits. In both cases, the energy of the 
leading parton is assumed to be large enough that its an- 
gular deflection can be neglected. One analytic limit ap- 
plies to thin plasmas where the mean number, n = L/X, 
of jet scatterings is small The other limit is the 

thick plasma one where n ^> 1. This mean number 

is a measure of the opacity or geometrical thickness of 
the medium: 
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For a constant (box) density p = N/(LA±), L is the tar- 
get thickness and A = l/(pa e i) is the average mean free 
path. For a sharp cylinder geometry we can interpret 
L = 1.2 A 1 / 3 = R s as the nuclear radius. For a more real- 
istic 3+1D expanding Gaussian cylinder with /j(xj_, t) = 
(T /T)p exp(-( X l + AT 2 )/R 2 G )l (2\x ± \AT/R 2 G ), L is 
replaced by the equivalent rms Gaussian transverse ra- 
dius Rq = 0.75 A 1 / 3 fm. The rightmost form in ([j]) is 
obtained by averaging over the expanding Gaussian cylin- 
der. Here To is the formation time of the plasma. 

At RHIC energies (y/s — 200 AGeV) the expected ra- 
pidity density of the gluons is dN/dy ~ 1000 for A = 200. 
With an elastic cross section <r e ; ~ 2 mb and a plasma for- 
mation time ~ 0.5 fm/c, we obtain h ~ 4. This suggests 
that neither analytic approaches may be strictly appli- 
cable in practice. However, because of the non-abelian 
analog jl|6) of the Landau-Pomeranchuk-Migdal (LPM) 
effect, the radiation intensity angular distribution and 
total energy loss are controlled by the combined effect 
of the number of scatterings n — L/X and a formation 
probability pf = L/lj = Lp 2 / (2xE) of the gluon in the 



medium. In the induced energy loss was shown to be 
nonlinear in the nuclear thickness, AE oc p 2 L 2 / X assum- 
ing the ti>l limit. In ||, we showed that the angular 
pattern was even more strongly nonlinear in the exclusive 
(tagged target) case. We show below that in the inclu- 
sive case, important features of the radiation pattern are 
in fact governed by the product npf = L 2 (p 2 /k + X). Be- 
cause of this, even the first order in opacity reproduces 
the L 2 dependence of the energy loss and the range of ap- 
plicability of the finite opacity expansion derived below 
may extend to realistic targets. 

Another important motivation for our approach is the 
apparent absence of jet quenching observed at SPS en- 
ergies ||[l0). As we find below, the opacity expansion 
shows that finite kinematic constraints suppress greatly 
the non-abelian energy loss at SPS energies. The work 
reported in this Letter extends Ref. [^) by including vir- 
tual corrections || necessary to compute inclusive (vs. 
tagged) jet quenching. Detailed derivation of these re- 
sults is given in Ref. [ pd| . 

Hard Gluon Distribution. At zeroth order in opacity, 
the gluon emission from the hard production vertex is 
approximately given by 
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where x = k + /E + sa oj/E, and Cr is the Casimir of 
the (spin 1/2) jet in the dn dimensional color represen- 
tation. For a spin 1 jet the gluon splitting function must 
be substituted above. The differential energy distribu- 
tion outside a cone defined by > p 2 is given by 
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where the upper kinematic limit is 

ki max = min [AE 2 x 2 AE 2 x{\ - x)] . (4) 
The energy loss outside the cone is then given by 
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in the leading \og(E/p) approximation (LLA). While this 
overestimates the radiative energy loss in the vacuum 
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(self-quenching) , it is important to note that AE^ ' / E ~ 
50% is typically much larger than the medium induced 
energy loss. 

Opacity Expansion. To compute the induced radiation, 
we assume as in [^,^) that the quark gluon plasma can 
be modeled by N well-separated, i.e. A 3> color- 
screened Yukawa potentials. In contrast to ||, we con- 
sider here the inclusive gluon distribution induced by a 
finite medium which remains unobserved. Without target 
tagging, we must add double Born (Virtual) amplitudes 
Vi to the real (Direct) amplitudes Dj derived in |5|. We 
denote by Go the basic eikonal hard emission amp itudc 
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where we use the shorthand notation of H, H = 
kj^/k^, ujq = k 2 L /2xE. The hard parton originates at zo, 
and c denotes the color matrix T c in its d R dimensional 
representation. In Ref. (ll) we developed an iterative 
procedure to generate from this amplitude the sum of all 
amplitudes to any order in opacity including both direct 
and virtual terms. 

To first order in opacity it is sufficient to consider 
the classes of amplitudes denoted Go£>i (sum of three 
graphs) and G0V1 (sum of four graphs) resulting from di- 
rect and virtual terms due to a scattering at position Z\ . 
The cross section for the induced radiation consists of 3 2 
real and 2x4 double Born contributions that sum to a 
simple result 

Tr (G Q D X D\G\ + (GoV^l + h.c.)) 

= 4g 2 C R C A d R [(-2d -Br) (1 - cos(cjiAzi)) ] , (7) 

with C m and uj m obtained from H and ujq through the 
substitution kj^ => k^ — q_i_ m and B m = H — C m . Note 
that the interaction of the radiated gluon in the medium 
brought in a factor Ca = N c . Unlike the tagged case [||, 
the above first order correction to the hard (factoriza- 
tion) distribution in Eq.(|^) has no simple classical cas- 
cade limit. In the tagged case, the result contained terms 
with color factors C R d R as well as C R CAd R that were 
easily interpretable in terms of jet rescattering and in- 
duced radiation. In the inclusive case, the virtual correc- 
tions lead to long range non-local interference effects that 
have no classical analog. On the other hand, a remark- 
able "color triviality" O] is found at the inclusive level 
where all jet rescattering effects cancel fl. As show in 
detail in Ref. |fl]], the color factor at order n — 2 is sim- 
ply C R C\d R , in contrast to the vastly more complicated 
color structure of the tagged case [|| . 

The non-abelian LPM effect is seen in Eq. (Q) as arising 
from the gluon formation factor 
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where Azi = Z\ — z . This must also be averaged over 
the longitudinal target profile, n(z). For a box den- 
sity § of thickness L = 1.2 A 1 / 3 , n = 9{L - z)/L. For 



analytic simplicity we take instead an exponential form 
n e (z) — e~ z / L " /L e . In order to compare results for the 
two cases, we must require identical mean target depths, 
i.e. L e = L/2. With n e (z) the ensemble averaged forma- 
tion factor is 
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The formation probability in this case is controlled by 
simple Lorentzian factors. 

Averaging over the momentum transfer qi_L via the 
color Yukawa potential leads finally to the gluon double 
differential distribution 



_ dN^ _ dN^ l_ r q ™ d2 

dxdk 2 ^ dxdk 2 ^ X g J 



H<£± + v 2 ) 2 

2k ± -q 1± (k-q 1 ) 2 L £ 2 
16x 2 £: 2 + (k-qi)|L 2 ' 



(10) 



where the opacity factor L/\ g — Na^ /A± arises from 
the sum over the N distinct targets. Note that the ra- 
diated gluon mean free path X g = {Ca/C r )\ appears 
rather than the jet mean free path. It is the color triv- 
iality of Eq. (0) that allows us to absorb C R factor into 
the hard distribution, Eq. (||), and Ca factor into \ g . 

The upper kinematic bound on the momentum trans- 
fer q 2 max = s/4 ~ 3£/i, (l/fi ff = l/^ 2 - 1/( M 2 + 9 2 nax )). 
For SPS and RHIC energies, this finite limit cannot be 
ignored as we show below. 

The second order contribution in opacity requires a 
more complex calculation involving the sum of 7 2 direct 
and 2 x 96 virtual terms as discussed in |ll| . The result 
is 

Tr ( GqD x D 2 d\d\gI + {G ViD 2 DlGl + h.c.) 

+ (GoD 1 V 2 d\gI + h.c.) + (GoViViGl + h.c. 

= Ag 2 C R C 2 A d R [2Ci-Bi (1 - cos(wiA«i)) 
+ 2 C 2 • B 2 (cos(w 2 Az 2 ) - cos(w 2 (Azi + Az 2 )) 

- 2 C( 12 ) • B 2 (cos(w 2 Az 2 ) - cos(cj (12 ) Azi + w 2 Az 2 )) 

- 2 C (12) • B 2(12 ) (1 - cos(w (12 )Azi)) ] , (11) 

where with C( m „) and u>t mv \ obtained from H and uiq 
through the substitution kj^ =>• kj^ — q_i_ m — q± ti and 
B m ( n n = C m — C(pi) ■ Note the color triviality of this in- 
clusive result j6| JT2j_ in contrast to the tagged target case 
in H . Both EqsT7|7|,|i"Tl) actually hold for either quark or 
gluon jets in the small x approximation. 

We note further that in the two important limits 
|kj_| — > and |kj_| — > oo the distributions vanish for 
any fixed x due to azimuthal angular averaging 
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This is clearly seen in Fig. 1. 

The average over the scattering points, Z\,z% in the 
second order case is performed with the exponential den- 
sity profile as follows: 



3 dA Zl f°° 3 dAz 2 



3(A Z1 +Az 2 ) 



(13) 



where Azi = — We must use L e = L/(n + 1) 

for n th order in opacity in order to insure that the first 
moment of the m th scattering center is identical to that 
in a box distribution {{z m ) n — mL/(n+ 1)) as discussed 
in §. 

Numerical results comparing the first and second or- 
der in opacity corrections to the hard distribution Eq. (j^) 
are illustrated in Figs. 1 and 2. We consider a 50 GeV 
quark jet in a medium with X g = 1 fm. A screening scale 
/i = 0.5 GeV and a s = 0.3 were assumed. The "angular" 
distribution in Fig. 1 shows that the first order angular 
distribution is wider than the medium independent hard 
(l/kj_) distribution. The second order corrections redis- 
tributes the gluons further. In Fig. 2 the integrated 
contributions to the gluon intensity, dl / dx, are shown for 
the same conditions as in Fig. 1 . The induced intensity is 
concentrated at small x in contrast to the relatively con- 
stant intensity originating from the hard "self-quenching" 
term (||). The long 1/x tail contributes however a loga- 
rithmic factor \og(E/fi) as we discuss further below. The 
second order correction suppresses the intensity at small 
x and enhances it somewhat at higher x > 0.2 . We note 
that these two effects tend to cancel in the integrated 
energy loss as seen in Fig. 3. 

Radiation Intensity and Energy loss. To gain analytic 
insight into the above numerical results, we consider the 
first order induced radiation intensity dE 1 ' /dx in the ap- 
proximation that k^_ max = oo. This allows us to change 
variables = k^ — qu_ in Eq. ( |l0| ) and express the inte- 
grand in the azimuthal <p integral as a partial derivative 
with respect to kj_. The remaining integral can be 
performed then analytically, resulting in 
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where 7 = L/1 2 / (4xE) is a formation probability and 
5 = fi/E is a measure of the importance of finite kine- 
matics. The formation function is given by 
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which in the 6 — > and 7 C 1 limit reduces to a simple 
form /(7, 0) w n-f oc L. It is this latter limit that leads 
to the characteristic quadratic dependence on L: 
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FIG. 1. The medium induced gloun differential multiplicity 
normalized by the factorization distribution eq.(Q) is plotted 
vs. k^/^ 2 for opacity L/\ g = 1, 2, 3. Solid curves show first 
order (|l^) and dashed including second order in opacity ( pd[ ) 
for a gluon with x = 0.2. (Quark jet of Ej et = 50 GeV and 
ft = 0.5 GeV). 
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FIG. 2. The contributions to the induced radiation inten- 
sity to first and second order in opacity are compared to the 
medium independent intensity (|^). The same conditions as 
in Fig. 1 are assumed. 
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This formula breaks down at both x — > and x — > 1 
because |kjj max cannot be approximated by 00 and be- 
cause the small x approximations used above break down 
as x — > 1. 

The total radiative energy loss is then given by 
C R a s L 2 [i 2 
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FIG. 3. The radiated energy loss of a quark jet with energy 
E jet = 5,50,500 GeV (at SPS, RHIC, LHC) is plotted as a 
function of the opacity L/\ g . (X g = 1 fm, /i = 0.5 GeV). Solid 
curves show first order, while dashed curves show results up 
to second order in opacity. The energy loss (solid triangles) 
from Ref. @ (with v = 2.5) is shown for comparison. 

with N(E) — 4 if the kinematic bounds are ignored as in 
Eq. (fl6l). In practice, we emphasize that finite kinematic 
constraints cause N(E) to deviate considerably from the 
asymptotic value 4. We find that N(E) = 7.3, 10.1, 24.4 
for E = 500, 50, 5 GeV. Together with the logarithmic 
dependence of energy, these kinematic effects suppress 
greatly the energy loss at lower (SPS) energies as seen in 
Fig. 3.. This is in contrast to the approximately energy 
independent result in Ref. Q where the finite kinematic 
bounds were neglected. 

Conclusions. We calculated the effect of final state 
interactions on the induced gluon differential distribu- 
tion up to second order in opacity for hard jets produced 
in nuclear reactions. This work generalizes Ref. ]5| by 
taking into account virtual corrections to calculate inclu- 
sive rates and provides an complementary analytic ap- 
proach to clarify nonlinear jet quenching effects predicted 
in Refs. §§§. The inclusive xdN/dxdk 2 ± , dl/dx, and 
AE, was studied as a function of nuclear thickness for 
jet energies in the SPS, RHIC and LHC range. One 
of our main results is the demonstration that the sec- 
ond order contribution ( pT| ) to the integrated energy loss 
remains surprisingly small up to realistic nuclear opaci- 
ties L/Xg ~ 5 ( except at low SPS energies). The rapid 
convergence of the opacity expansion even for realistic 
opacities results from the fact that the effective expan- 
sion parameter is actually the product of the opacity and 
the gluon formation probability L/j, 2 jxE (see (|l5|)). The 
leading quadratic dependence of the energy loss on nu- 
clear thickness therefore arises from the simple first order 
term (|7|) in this approach. The detailed pattern of angu- 
lar broadening and the x dependence is also dominated 
by the first order contribution. We note that our first 
order (effectively power-law) |kj_| and x distributions, 



however, differ considerably from the Gaussian form ob- 
tained in the eikonal resummation approach of |]|| that 
applies to thick targets. Surprisingly, on the other hand, 
the magnitude of the integrated AE(V) is rather similar 
at RHIC energies. 

At SPS energies kinematic effects suppress greatly the 
energy loss relative to || . Our estimates provide a natu- 
ral explanation for the absence of jet quenching in Pb+Pb 
at 160 AGeV that has been a puzzle up to now [p|,^0[. 
The short duration of the dense phase further limits the 
effective opacity at the SPS. A duration of L/X g ~ 2 for 
example leads to a total energy loss only ~ 100 MeV, 
which is much too small to be observable in soft mul- 
tiple scattering background pfj[ . At RHIC energies, on 
the other hand, a significant nonlinear (in A) pattern of 
suppression [Q,^| of high p± hadrons relative to scaled pp 
data should be observable to enable a direct test of non- 
abelian energy loss mechanisms in dense matter. We note 
finally that the simplicity of the first and second order re- 
sults (0,|ll|) will make it possible to improve significantly 
Monte Carlo event simulations of jet quenching 
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[Note added in proof: U. Wiedemann informed us of 
an opacity expansion derived independently in p3|.] 
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